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Abstract 

We have studied the specific heat of the {Ns + Nb) model for an Ns-hody harmonic oscillator 
(HO) system which is strongly coupled to an A'^^-body HO bath without dissipation. The system 
specific heat of Cs(T) becomes Nsks at T — t- oo and vanishes at T = in accordance with the 
third law of thermodynamics. The calculated Cs(T) at low temperatures is not proportional to 
Ns and shows an anomalous temperature dependence, strongly depending on Ns, Nb and the 
system-bath coupling. In particular at very low (but finite) temperatures, it may become negative 
for a strong system-bath coupling, which is in contrast with non-negative specific heat of an HO 
system with Ns = 1 reported by G-L. Ingold, P. Hanggi and P. Talkner [Phys. Rev. E 79, 061105 
(2005)] . Our calculation indicates an importance of taking account of finite Ns in studying open 
quantum systems which may include an arbitrary number of particles in general. 

PACS numbers: 05.30.-d, 05.70.-a, 65.40.Ba, 65.80.+n 



*hideohasegawa@goo.jp 



1 



I. INTRODUCTION 



In recent years, there has been considerable interest in physical properties of small-scale 
systems, which are prepared by advanced new techniques pj. Stimulated by this devel- 
opment, a study on small-scale systems which is one of important areas in classical and 



quantum statistics, has been extensively made Theoretical studies on open systems 
have been made with the use of the Caldeira-Leggett (CL) type models 0-7], in which a 
single particle is assumed to be subjected to a bath consisting of uncoupled harmonic oscilla- 
tors (HOs). The CL model was originally proposed for an infinite bath {Nb — )■ oo). Recent 



studies, however, have employed the C 
small system coupled to finite bath 



model with finite A''^ for studies of properties of 
hJ. Refs. 0,3 have studied thermalization of a 



particle (the system) coupled to a finite bath, showing that a complete thermalization of the 
particle requires some conditions for relative ranges of oscillating frequencies in the system 
and bath. The energy exchange between particles in a rachet potential (the system) and 
finite bath (A^^ = 1 — 500) has been investigated [lo|. 

In CL-type models having been proposed for open systems , there are two issues 
when they are applied to realistic open systems: (a) the number of particles in a system is 
taken to be unity {Ns = 1) and (b) a system-bath coupling is assumed to be weak although 
the overall damping can be strong. As for the issue (a), a number of particles in a system 
is required to be finite since a generic open system may contain any number of particles. 



CL-type models with Ns = 2 have been investigated in Refs. [Uj, ll2|. In our previous 



paper 



13| . we proposed the (A^^ + A''^) model in which a finite A^5-body system {Ns > 1) is 



coupled to finite A'^g-body bath. It has been shown that calculated energy distributions of 
a system show intrigue properties as functions of Ns, Nb and a system-bath coupling 
As for the issue (b), Refs. 



3. Q, Q have pointed out ambiguities in defining physical 



quantities such as energy and specific heat when a system-bath coupling is not weak. Two 
different routes toward the evaluating of a system energy have been proposed for a system- 
plus-bath Q, [isj. The first route is based on the system partition function as given by 
a ratio between the total and bath partition functions [Eg. ( I45p ]. which is traditionally 
identified as the partition function of an open system 16|. The second route is based on 
the expectation value of the system Hamiltonian averaged over the total Hamiltonian. In 
the limit of vanishing or weak interaction, the two definitions yield the same results. It 
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is, however, not the case in general when the interaction is not weak (finite). The specific 
heat of the system of a single free particle coupled to a bath described by the Drude model 
has been studied with the use of the two evaluation methods [l4, [isl, 17|. Specific heats 



obtained by the two routes are different not only at low temperatures but also in the leading 
high-temperature correction terms Q, [isl . In particular, the specific heat in the first route 
is negative at low temperature while that in the second route is not. In contrast, similar 
calculations for a single HO system subjected to a single bath oscillator or Drude bath 
have reported that the first route does not yield a negative specific heat [l^. The difference 
between negative specific heat in a free particle system and non-negative one in an HO system 
within the first route is attributed to the fact that the degree of freedom in the former is 
smaller than that of the latter [17]. The obtained negative specific heat is not unphysical 
because the system specific heat calculated by the first route should be interpreted as a 
change in the specific heat of the environment when a system degree of freedom is attached 
0, Q, We should note that studies of Refs. flQ Q have been made for the CL 
model with A^^ = 1. It is not clear whether conclusions obtained in Refs. Iisl . are 
valid in a general case, for example, for finite Ns (> 1). 

It is interesting to study thermodynamical properties of the {Ns + Nb) model proposed 
in Ref. [isj, which is worthwhile for us to get some insight to the issues (a) and (b) and 
which is the purpose of the present paper. In conventional studies, deterministic dynamics 
of particles in a system-plus-bath is replaced by the stochastic Langevin equation. By using 
the alternative method in this study, we will evaluate eigenfrequencies of the {Ns + Nb) 
model where the system is not dissipative for finite Nb 18[]. The energy and specific heat 
have been calculated with the use of the first route for evaluating energy mentioned above. 

The paper is organized as follows. In Sec. II, we briefi y ex plain the {Ns + Nb) model 
for coupled HO system subjected to uncoupled HO bath [l3|, emphasizing a difficulty in 
solving the A^^-coupled quantum Langevin equations. With the use of canonical transfor- 
mation, we derive the feasible quantum Langevin equation, with which the expression for 
eigenfrequencies of the model is obtained. Eigenfrequencies are analytically evaluated for a 
bath described by the identical-frequency model {oOn = ooq ioi n = 1 to A''^) Thermo- 
dynamical quantities are expressed in terms of obtained eigenfrequencies. We have made 
numerical calculations of temperature dependences of the energy and specific heat, changing 
parameters of Ns, Nb and a system-bath coupling. An origin of the negative specific heat at 
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low temperatures is illustrated. In Sec. Ill we study similar two models for uncoupled HO 
system subjected to uncoupled HO bath (Sec. IIIA) and for coupled HO system subjected 
to coupled HO bath (Sec. IIIB). The final Sec. IV is devoted to our conclusion. 



II. ADOPTED (Ns + Nb) MODEL 

A. Coupled quantum Langevin equations 

We consider the {Ns + Nb) model in which the a one-dimensional A^^-body system (Hs) 



is subjected to an NB-hodj bath {Hb) by the interaction [Hj) 13|]. The total Hamiltonian 
is assumed to be given by 



H = Hs + Hb + Hj, (1) 



with 



'Pl DQl K.^ ^ , 



k=l 



(2) 



= + (3) 

n=l ^ 

Ns Nb 

Hi = EEl^(^^-^-)'' (4) 

fc=l n=l 

where Pk (pn) and Qk (qn) express the momentum and position operators, respectively, of an 
HO with mass M (m) in the system (bath), D and K denote force constants in the system, 
and Ckn is a system-bath coupling. The system is subjected to a bath given by Eq. 
consisting of a collection of uncoupled HOs with oscillator frequencies of {un}- Operators 
satisfy commutation relations, 

[Qk,Pe] = ifiSke, [qn,Pm] = ih6nm, [Qk,Qe] = [Pk,P£] = [qn,qm] = [Pn,Pm] = 0. (5) 

In a general case of -D 7^ and K 0, Eq. ([2]) expresses the system including coupled 
HOs. In the limit of 7^ and K = 0, the system consists of a collection of uncoupled 
(independent) HOs. 

In conventional approaches to the system-plus-bath model, we derive the quantum 
Langevin equation, employing the Heisenberg equation, 

ihO = [0,H], (6) 



where O stands for an arbitrary operator and a dot denotes a derivative with respect to 
time. Equations of motion for Qk and g„ are given by 

Nb 



MQk = -DQk - K{2Qk - Qk-i - Qk+i) ~^Ckn{Qk ~ Qn), (7) 

n=l 

Ns 

mqn = -mulqn + ^ CknQk, (8) 



fc=i 



with 



Ns 



It 11 



+ Ckn. (9) 

Substituting a formal solution of g„(t) into Eq. ([7]), we obtain the quantum Langevin 
equations given by jisl 

Ns 

MQk{t) = -DQk{t) - K [2Qk{t) - Qk-i{t) - Qk+i{t)] - MY^ikiQiit) 

1=1 

Ns j-t Ns 

-J2 ^^'^^ - t')Q^it') dt' -J2^MQi{0) + (kit) for A; = 1 to Ns, (10) 

with 

iVfl / N 

M^ki = y](ckjki-'-^], (11) 

lkl{t) = I I COSWnt, (12) 

Nb / . \ 

Cfc(^) = V'cfcn I g„(0) cosa)„t + — sinwnt). (13) 
n=i V a;n / 

Here ^ke denotes the additional interaction between k and ^th particles in the system induced 
by couplings {cfc„}, 'jkiit) stands for the memory kernel and (k is the stochastic force. By 
using averages in initial values of qn{0) and gn(0), 

{m^lqniOr)B = m{UOr)B = (^) coth , (14) 

we obtain the fluctuation-dissipation relation, 

lmMt') + Q{t')Ut))n = (^)coth(^)cosa;.(t-t'), (15) 

^ kBT^kiit - 1') for /3 -> 0, (16) 
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where {■)b signifies the average over initial states of the bath. In the case of Ns = 1, $,kk in 



Eq. f[TT|) expressing a shift of oscillator frequency due to an introduced coupling, vanishes 
if we adopt c„ = mw^ [4, 5|. In the case of A^^ 7^ 1, however, it is impossible to choose 
{ckn} such as C,ke = for all pairs of {k,i), then Q^. is inevitably coupled with Qi {i 7^ k). 
Because of these couplings between HOs, the A^5-body system cannot be simply regarded 
as a sum of systems with Ns = 1. Although the quantum Langevin equations given by 
Eqs. f|T0|) - f|T3|) are formally exact, it is difficult to solve them because they are given by 
A^5-coupled integrodifferential equations. 



B. Quantum Langevin equation with canonical transformation 

It is possible to derive the quantum Langevin equation which has a simpler structure 
than that given by Eq. (|T0 l) -( |T3l) . We assume that Ns is even without a loss of generality. 
Imposing a periodic boundary condition. 



Q 



Ns+k 



Qki PNs+k — Pki 



we employ the canonical transformation [2 



7Vs/2-l 



Qk 



= E 



/Ns 



^i{2TTks/Ns)Q^^ 



(17) 



(18) 



-Ns/2 
Ns/2-l 



s=-Ns/2 

Note that the boundary condition is satisfied in Eqs. ( |T8l) and ( 1T9|) and that the set 
1(1/ \/Ns) e^'-^^^/^'^^^ I is orthogonal and complete in a periodic domain of the oscillator 



label k 
becomes 



20| . With the canonical transformation given by Eqs. (fT8|) and ( !T9|) . Hs in Eq. 



Ns /2-I 

E 

s=-Ns/2 



p:Ps , MnlQiQ, 



2M 



+ 



with 



ns 



D + AKsin^ ( ^ ) for . = + 1, .., ^ - 1, 



(20) 



(21) 



where the commutation relations: 



(22) 



hold with Q* = Q_, and P; = P_,. 

The canonical transformation given by Eqs.f llSp and (|T9|) with an assumption, 



(23) 



leads to Hi given by 



TVs A^s/2-1 ^ ^ ^ 



(24) 



n=l 



n=l 



n=l 



From Eqs. ([3]), f l20|) and flMj) . we may derive equations of motion for Qs and 



m4o 



for s 7^ 0, 



-Mli^Qo + ^ c„g„ for s = 0, 



71=1 



mg„ = -mujlqn + a/ NsCnQ. 



with 



Nb 

Mnl + Cn for s 

n=l 



Ns 
2 ' 



TVs 
2 



' 2 ' 



mw^ = mul + NsCn, 



(25) 
(26) 
(27) 

(28) 
(29) 



which show that effective frequencies of Clg and Ur are increased with increasing the system- 
bath coupling c„. 

Substituting a formal solution of qn(t) of Eq. ( 1271) into Eq. ( l26i) . we obtain 



MQo(t) 



where 



Mfi2_^(0) go(t)- / i{t-t')Qo{t')dt' + C{t)-^{t)Qo{Q), (30) 





(31) 
(32) 



7(t-t') = NsJ2 I COSUr.it -t') 
n=l \"^^n 

Nb 



at) 



n=l 



g„(U) cos LOnt H z — smw„t 



Equations (I30l)-(l32l) express the desired quantum Langevin equation for Qo, an s = com- 
ponent of Qs, which has simpler structure than those given by Eq. (fT0l)-(IT3l). Note that 
solutions of Qs with s 7^ are obtained from Eq. (!25|) . Once Qs for all s are obtained, Qk 
is obtainable with the use of the canonical transformation given by Eq. (fT8|) . 



We note that the quantum Langevin equation given by Eq. (l30ll is similar to that in 
the conventional CL-type model with Ns = 1 where g„ couples to a single particle in 
the system. We may employ various methods such as the quantum Langevin and master 
equations, which have been widely adopted for the CL-type model. In the case of Nb — oo, 
a sum in the kernel given by Eq. fl31l) is converted into the integral and the kernel 7(t) may 
be approximated by the Drude or Ohmic one. 



C. Eigenfrequencies of small open systems 



In the present study, we pay our attention to small systems subjected to non-dissipative 
finite bath ^ ^ 

ft 



181 ]. Equation (l3Ul) - fl5^ may be rewritten as 

MQoit) = -MQlQo{t)+ [ fi{t-t')Qo{t')dt' + C{t), 







with 



Mt-t')^iV.X:(-|-)sm^.(t-t')^^7(t-t'), 



(33) 



(34) 



which are more tractable than Eqs. (!30l) -( j32l) for finite non-dissipative baths. Applying the 
Laplace transformation to Eqs. ( I33l) and ( l34l) . 



e-''Qo{t) dt, 
e'^^Qnit) dt, 



we obtain 



Qo{z) = Giz) 



NsCn /g„(0) + 2;g„(0) 



with 



G{z) 



n=l 



M 



Z"^ + LOI 



(35) 
(36) 

(37) 



(38) 



^ mM{z'^ + Col) 

We note in Eqs. ( I25l) -( l271) that Qs for s 7^ are decoupled from the rest of variables, 
while Qq is coupled with g„. (A^^ — 1) eigenfrequencies of Qsif) are given by VLs in Eq. (125|) . 
Remaining [Nb + 1) eigenfrequencies are given by roots of G^—w)'^ = 0, 



V 



VlI-NsY. 



^ ■mM{y'^ — cj^) 



0. 



(39) 
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Alternatively, we may obtain {Nb + 1) eigenfrequencies as follows: Calculating the de- 
terminant derived from Eqs. fl26|) and f l27|) which is expressed in a basis of (s = 0, n = 
1,2,3,...), 



-C2VNs -csVNs ... 
-ciViVs m(w2-z/2) 
-caViVs miCj^-u^) 



0, 



and using a sweeping method to make a triangle determinant, we obtain an equation for 
eigenfrequencies , 



n=l 



miujl — p"^) 



n=l 



(40) 



which is equivalent to Eq. (15^ . 

It is difficult to analytically solve Eq. ( 139|) or (HOj) in a general case of {wn}, which requires 
numerical methods. However, when we adopt the identical-frequency model for the bath 



Un = cuo, c„ = c for = 1 to A''^, 
we may easily obtain eigenfrequencies of z/j (z = 1 to Ns + A'^), 



(41) 



i 


1 




Nsl2 + 1 






A^s + 1 




ATcj + NbII + 1 








"-Afs/2 








"Afs/2-l 


^0 











with 



~,2 







MVlI + A^B c, 
+ Ns c. 



(fig - ^2)2 + 



4A^sArgc2 
Mm 



(42) 

(43) 
(44) 



where fig is given by Eq. (12T]) . It is easy to see that in the limit of no couplings (c = 0), 
eigenfrequencies are given by 



% 


1 




Nsll + 1 




A^s 


A^5 + l 




Ns + ArB/2 + 1 




A^s + iVfi 




"-TVs/a 








"Afs/2-1 






^0 




u;2 
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In the case of Ns = Nb = 1 with K = 0, D = fs, uq = and c = fs, Eq. (|l2l) yields 
eigenfrequencies which are equivalent with those expressed by Eq. (17) of Ref. [l7|, related 
discussion being given in Appendix A. In a pedagogical case of Ns = 2 and Nb = 1, we may 
exactly solve a system-plus-bath, transforming original variables of Qk and Pk {k = 1,2) 
into center-of-mass and relative variables. Obtained eigenfrequencies agree with those given 
by Eq. (H2|) as is shown in Appendix A. 



D. System partition function, energy and specific heat 

It has been shown in the preceding subsection that the system-plus-bath may be effec- 
tively expressed as a collection of independent HOs with eigenfrequencies {z^i}. The system 
partition function is expressed in terms of eigenfrequencies by 



with 



Ns+Nb 

Tre^^= n 

i=l 
Nb 



(45) 



Zb = Ttb e-^^- = n 



2sinh(/3nz/i/2)_ 
1 

2smh{Phuj/2)_ 



(46) 
(47) 



where Tr and Tr b denote a full trace over all variables and a partial trace over bath variables, 
respectively, and a bath frequency is given by Uj = Uq. The energy and specific heat of the 
system are given by 

dlnZs 



E, 



dp ' 



i=l 
Ns+Nb r 

1=1 
C-Cb 



, coth ( ^ 
2 V 2 



coth 



E 



phui 



2smh{(3hui/2) 



-kBNB 



2siTi\i{(3nujQ/2) 



(48) 
(49) 

(50) 
(51) 



where the first term of C expresses the specific heat of a system-plus-bath and the second 
term oi Cb denotes the specific heat of the bath alone. In the zero- and high-temperature 
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limits, Es and Cs become 



Es 



Cs 









Nsks for ksT ^ hui, Hwq. 



for ksT = 



for ksT ^ hui, hw^ 



for /c^T = 0, 



(53) 



(52) 



E. Model calculations for coupled HO systems subjected to uncoupled HO baths 

We have made numerical calculations of Es and Cs, by using Eqs. ( 149|) and (l50l) . The 
system-bath coupling c is assumed to be given by jlj] 



such that the interaction term in Eq. (jl]) including summations over Ylk=i Sn=i yi^ld 
finite contributions even in the limits of A^^ — ?■ oo and/or Nb — oo. Our (A^^ + A''^) model 
includes eight parameters: M, D, K, Ns, m, uq, Nb and cq. Bearing in mind coupled HO 
system subjected to uncoupled HO bath, we have employed K = 1.0, D = 0.0, m = 1.0, 
M = 1.0 and Uq = 1.0 otherwise noticed, related discussion being given in Sec. 111. A unit 
of energy of the model is given by eo = hwo- 

1. Isolated system and hath 

First we show calculated results of isolated system and bath with cq = 0.0. Figure [U^a) 
shows Eb/Nb and Es/Ns for Ns = 4, 10, 100 and 1000. The bath energy is Eb/Nb = 
at T = and gradually approaches /c^T at high temperatures. The Ns dependence of the 
system energy Es/Ns is not so evident except for very low temperature (below). Figure 
[U^b) shows CB/kBNB and Cs/kBNs for various A^^. The bath specific heat CB/kBNB 
follows Einstein's formula showing exponential decrease at low temperatures. In contrast, 
the system specific heat becomes Cs/kBNs = l/Ns at T ~ 0.0, although it reduces to zero 
at T = as shown by Eq. (I53k). This transition is more clearly seen when we introduce 
an infinitesimal coupling (Fig. [2]). At high temperatures, Cs/kBNs approaches unity as 
expected. In the intermediate temperature range, Cs nearly follows the linear T. The 
linear-T specific heat is examined in the inset of Fig. Wih) where Cs/kBTNs is plotted as 




c = 



NsNb 



(54) 
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a function of T. We note that Cs/ksTHs for Ns = 100 and 1000 is almost constant at 
0.01 < T < 0.6. It is noted that the hnear-T behavior of the calculated specific heat is a 
consequence of the one-dimensional coupled HO model adopted in this study. If we employ 
the three-dimensional coupled HO model for the system, we obtain the T^-specific heat at 
low temperatures. 

2. Effects of Co 

Figure [2] shows Cs/ksNs for various coupling strengths with A^^ = 4 and A^^ = 100. The 
system specific heat for cq = 0.0 is Cs/ksNs = 1/Ns at T ^ 0.0 as mentioned above. When 
a small interaction is introduced, Cs clearly reduces to zero at T = 0.0. With increasing 
the coupling strength, Cs is furthermore decreased at low temperatures where it shows an 
anomalous temperature dependence. For cq > 2.0, the specific heat becomes negative, and 
magnitudes of negative dips are increased with increasing cq. 

3. Effects of Ns 

The Ns dependent specific heat is shown in Figs. El^a), (b) and (c) where Cs/ksNs 
is plotted for cq = 0.1, 1.0 and 5.0, respectively, for various Ns with a fixed Nb = 100. 
Figure [3](a) shows that for cq = 0.1 an anomalous bump in Cs at low temperatures is 
gradually decreased with increasing Ns and it well follows the linear-T behavior at higher 
temperatures. For cq = 1.0 the temperature dependence of the specific heat is almost 
independent of Ns as shown by Fig. El^b). In contrast. Fig. El^c) shows that Cs strongly 
depends on Ns for cq = 5.0, for which magnitude of negative specific heat is much increased 
for smaller Ns- 

4. Effects of Nb 

The Nb dependence of the system specific heat is shown in Fig. HI where Cs/ksNs is 
plotted for cq = 0.1, 1.0 and 5.0 with various Nb and a fixed Ns = 4. In the case of cq = 0.1, 
the calculated specific heats are almost independent of A''^. In the case of Cq = 1.0, a bump 
in Cs is gradually decreased with increasing Nb- On the contrary in the case of cq = 5.0, 
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magnitudes of negative dips in the specific lieat are more significant with increasing Nb 
although the result for Nb = 1000 is nearly the same as that for Nb = 100. 



5. Origin of the negative system specific heat 

We will elucidate the physical origin of the negative system specific heat for a typical 
case of Ns = 4, Nb = 10 and Cq = 5.0, whose result has been presented in Fig. HI Chain 
and dashed curves in Fig. [S] express a total specific heat C and a bath contribution Cb, 
respectively, which arise from the first and second terms in Eq. flSUj) . The solid curve denotes 
the system specific heat of Cs {= C — Cb), which is given by the difference between the 
chain and dashed curves. The inset of Fig. [5]shows eigenfrequencies {coj} for cq = 0.0 (open 
circles) and Cq = 5.0 (filled circles). In Eq. fl50|) C is expressed in terms of {uj} for cq = 5.0 
which become larger than those for cq = 0.0 by an introduced system-bath interaction. 
As a consequence, the specific heat of system-plus-bath is suppressed at low temperatures 
compared to Cb for the bath HOs which is expressed in terms of uq. Then the system 
specific heat Cs given by the difference of C — Cb becomes negative at low temperatures. 
The obtained negative specific heat is not related with an instability of the system. 



III. DISCUSSION 



A. Uncoupled HO systems subjected to uncoupled HO baths 



In the preceding section, we considered coupled HO system with D = 0.0 and K = 1.0. 
On the other hand, when we adopt D = 1.0 and K = 0.0, the model given by Eqs. 
expresses the uncoupled HO system subjected to uncoupled HO bath. Eigenfrequencies of 

Ui {i = 1 to Ns + Nb) are given by 



i 


1 




Nsll + 1 




Ns 


Ns+1 




Ns + Nb/2 + 1 




Ns + Nb 




"0 




-\ 














^2 



with 



\ ^l+^l±\im-^l)' + 

D + Nbc, 
uiUq + Ns c. 



ANsNbc^ 
Mm 



(55) 

(56) 
(57) 
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In the limit of no couplings (c = 0), eigenfrequencies are given by 



i 


1 




Nsll + 1 










A^s + + 1 




Ns + Nb 




"o 




"o 






^0 








^0 



The system specific heat calculated with the use of Eq. (1501) is shown in Fig. El^a) where 
Cs/ksNs is plotted for various Ns with A''^ = 100 and cq = 5.0. For cq = 0.0, Cs/ksNs 
follows the Einstein specific heat [see the dashed curve in Fig. 1(b)]. When cq is much 
increased, the specific heat becomes negative at low temperatures. It is noted that the solid 
curve for Ns = 1 and A''^ = 100 shows a negative C5 at < ksT/eo < 0.805 (eo = ^0), 
while Cs for A^^ = Nb = 1 is non- negative 17|. Figure 6(a) clearly shows 



Cs{T;Ns,Nb) ^ NsCsiT;l,NB] 
= NsCs{T;1,Nb] 



for T ~ 0.0, 

for T — )■ 00, 



(58) 
(59) 



where Cs(T; Ns, Nb) denotes the system specific heat of the (Ns+Nb) model at temperature 
T. 

In order to examine an origin of the negative specific heat, we plot C, Cb and Cs for 
Ns = 4, A''^ = 10 and Cq = 5.0 in Fig. MJo), whose inset shows eigenfrequencies for cq = 0.0 
(open circles) and cq = 5.0 (filled circles). We obtain a negative at < kBT/e^ < 0.377 
where C < Cb although C5 > at kBT/e^ > 0.377. Eigenfrequencies for cq = 0.0 are 
Ui = 1.0 for z = 1 to 14, and they are increased by an introduced system-bath coupling of 
Co = 5.0 as shown in the inset. From a comparison between Figs. 5 and 6, we note that a 
negative specific heat is realized both in coupled and uncoupled HO systems subjected to 
uncoupled HO baths, independently of eigenfrequencies of the system. 



B. Coupled HO systems subjected to coupled HO baths 



We have so far considered that a system is subjected to a bath including a collection of 
uncoupled HOs. Here we will study a case in which a bath consists of coupled HOs. The 
system- plus-bath is described by the Hamiltonian given by Eqs. ([I])-(11D but Hb is replaced 
by 



Hp 



Nb 

E 

n=l 



— + —(a 
2m + 2 



qn+i) 



(60) 
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where kn stands for a force constant between neighboring particles in the bath. We assume 
that Nb is even and kn = k for = 1 to Nb, imposing the periodic boundary condition 
given by 



QNB+n — Qn^ PNB+n — Pn- 



(61) 



By using the canonical transformation with the identical-frequency model [Eq. (HT]) ]. we 
obtain eigenfrequencies given by (detail being given in Appendix B) 



i 


1 




iV5/2 + 1 




Ns 


Ns + l 




Ns + Nb/2 + 1 




Ns + Nb 








-I 




"TVs/a-i 






Z/2 







where 



with 



ANsNbc^ 
Mm 



(62) 



D + Mnl + cNb, 
muol + cNs-, 



AK I Txs 
M \Ns 

4k ( vrr 

sm 

m X^B 



= \ sin — 



for s 



Ns Ns I 1 Ns 1 

2 ' 2 ' '2 ' 




for r = + 1 ■ ■ ■ ^ - 1 



(63) 
(64) 

(65) 
(66) 



Effective frequencies of Qg and are increased with increasing c. In the case of no couplings 
(c = 0), eigenfrequencies are given 



i 


1 




Ns/2 + 1 




Ns 


iV5 + l 




+ + 1 




A^5 + iVs 




^-Ns/2 




"o 




02 

"Afs/2-1 


^-iVs/2 




^0 




^ !2 

^iVfl/2-1 



The system specific heat is expressed in terms of obtained eigenfrequencies {t'i}, 



Cs 



Ns+Nb 
i=l 

C — Cb, 



2sinh(/3nz/i/2) 



Nb /2-l 
r=-NB/2 



2smh{l3hur/2) 



(67) 
(68) 



where C and Cb express the first and second terms, respectively, of Eq. ( 1671) . System 
specific heats in zero- and high-temperature limits are given by Eq. (l53l) . 
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We have calculated the system specific heat with k = 1.0, K = 1.0 and D = 0.0, for 
which a unit of the energy is given by eo = Hy^k/m. Figure [Tl^a) shows the temperature 
dependence of Cs/ksNs for various Ns with A''^ = 100 and cq = 5.0. For cq = 0.0, Cs is 
positive, reducing to zero at T = 0.0 [see Fig. 1(b)]. When cq is introduced, C5 at T ~ 0.0 
becomes negative. Magnitudes of negative Cs are increased with decreasing Ns- In order 
to study an origin of the negative Cs, we plot C, Cb and Cs for A^^ = 4, A''^ = 20 and 
Co = 5.0 in Fig. [7|(b), whose inset shows eigenfrequencies for cq = 0.0 (open circles) and 
Co = 5.0 (filled circles). For cq = 0.0, eigenfrequencies of bath and system display dispersion 
relations given by Eqs. ( l65l) and ( l66i) . By an introduced system-bath coupling of cq = 5.0, 
eigenfrequencies are increased as shown by filled circles in the inset. This modification in 
eigenfrequencies leads to a negative C5 at < /c^r/eo < 0.386 where C < Cb- 



IV. CONCLUSION 



The specific heat has been studied of a small quantum system consisting of A^^-body HOs 
which is strongly coupled to an A'^g-body non-dissipative HO bath with identical frequency 
[Eq. fHT]) ] or with dispersed frequencies [Eq. f l66|) ]. The obtained results are summarized as 
follows: 

(i) Although the system specific heat of Cs{T) is proportional to Ns (extensive) in the 
high-temperature limit, it is not (non-extensive) at low temperatures except for a vanishing 
system-bath coupling, and 

(ii) Cs(T) displays an anomalous temperature dependence at low temperatures where it 
may become negative for a strong system-bath coupling. 

The item (i) implies that it is necessary to take into account finite A^^ for a study of small 
quantum open system, although previous studies have been made by exclusively using the 
CL-type models with A^^ = 1. The item (ii) is in contrast to the result of Ingold, Hanggi 
and Talkner p^] who reported that the specific heat of an HO system with Ns = 1 cannot 
be negative although that of a free damped particle system may be negative. They claimed 
;hat it is due to a larger specific heat of an HO by a factor of two than that of a free particle 
17 1 . The elucidation of non- negative specific heat of a HO system in Ref. [l7| which is valid 
for Ns = 1, cannot be applied to the case of arbitrary Ns- 

The system specific heat of Cs = C — Cb given by ( 1^ or (jHE]) shows that Cs expresses 
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a change of the specific heat when the heat bath is enlarged by couphng it to system degree 
of freedom. The appearance of a negative specific heat in the item (ii) is attributed to 
increased eigenfrequencies when the system is coupled to the bath, by which the specific 
heat of C is suppressed as C < Cb and then < at low temperatures. The expression of 
Cs = C — Cb reconciles with the experimental procedure to determine the specific heat of 
the system attached to the bath. First one measures the specific heat of the empty container 
and then subtracts this value from the measured specific heat of the combined system-plus- 
bath to finally obtain the specific heat of the system |l4, 1^]. Specific heat anomalies 
which arise from an attachment of the system to the environment may be common in small 
open quantum systems. It would be interesting to examine the items (i) and (ii) by relevant 
experiments for small-scale systems. 
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Appendix: A. Eigenfrequencies of (1 + 1) and (2 + 1) models 



1. The (1 + 1) model 



From Eqs. ([I])-(!4j), the Hamiltonian for the case of Ns = Nb = 1 is given by 



H 



2M 



+ 



2 ^ 2m ^ 



— + 5W-^)^ 



for which we easily obtain eigenfrequencies. 



D + c muQ + c\ I f D + c muQ + 4c^ 



M 



+ 



m 



± 



M 



m 



+ 



Mm 



The CL-type model employed in Ref. 17| is given by 

)2 /\/rn2/o2 p2 
2m 



H' 



f fB, 



2M 2 2m 2 

which is related with the Hamiltonian of Eq. flAl|) by Z) = Mfi^, c 
Equation (1A2p is equivalent to Eq. (17) in Ref. jl7|. 



(Al) 



for i = 1,2. (A2) 



(A3) 

fB and Uq = 0. 
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2. The (2 + 1) model 



The Hamiltonian for the case of Ns = 2 and Nb = 1 is given by H = Hs + Hb + Hj with 

Hs = ^{Pl + P|) + y(g? + Ql) + ^(Qi - Q2)^ (A4) 

HB = ir.+ (A5) 

(A6) 



i^/ = ^[{q-Qif + {q-Q2f], 
meanings of all terms being trivial. Introducing center-of-mass and relative variables, 

Qc = ^iQl+Q2), Qr = Ql-Q2, 



we obtain Hs and Hj given by 

Hs = 



P1+P2, Pr = -{Pl-P2), 



PI 2 Pt DQl 



where M^. = 2M and = M/2. Equations of motion for Qc, Qr and q are given by 

McQc = -2{D + c)Qc + 2cq, 

MrQr = --{D + 2J + c)Qr, 

mq = —{muQ + 2c)q + 2cQc- 



(AT) 
(A8) 

(A9) 
(AlO) 

(All) 
(A12) 
(A13) 



We note that Qc is coupled with q whereas Qr is decoupled from Qc and q. A simple 
calculation leads to three eigenfrequencies given by 



±{D + 2J + c) 



±d \— I + 



for i = 1, 
for z = 2, 3. 



(AM) 



Equation f lA14p agrees with Eq. (jl2]) for A^5 = 2 and Nb = I because Mfig = D + c, 
itlCoq = itluq + 2c and J = 2K which is due to a double counting of interactions for Ns = 2: 

K ELiiQk - Qk+if = K[{Q^ - Q^f + {Q2 - Qif] = J{Qi - Q2f. 

It is evident from Eqs. ([50]), fjA2l) and flA14D that we obtain 



C5(T; 2,1) ^ 2C5(T;1,1) for T ~ 0.0, 

= 2 Cs{T] 1, 1) = 2kB for T ^ cx). 



(A15) 
(A16) 
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where Cs(T; 1, 1) and Cs(T; 2, 1) denote T-dependent system specific heats for the (1 + 1) 
model and (2 + 1) model (with J = 0.0), respectively. 



Appendix: B. Baths consisting of coupled harmonic oscillators 

We consider an A^^-body coupled system subjected to an NB-hodj coupled bath, whose 
Hamiltonian is given by Eqs. dl]), (E]), (H]) and flBU]) . Applying the canonical transformation 
given by Eqs. ( ITS]) and ( TT^ to the system with the boundary condition given by Eq. (IT7|1 . 
we obtain Hs given by Eqs. fl2U]) and fl^T]) . When we apply the canonical transformation 
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20|, 



r=-NB/2 
Nb/2-1 

p^= J2 e^^2-"^/^-V~„ (B2) 



^ r=-iVfl/2 



to the bath with the periodic condition given by Eq. fl6T]) . Hb in Eq. fl60|) becomes 



Nb/2-1 / 2~*- \ 

r=-NB/2 ^ ^ 



with 



4A;\ 



"'=iij="H£j forr = -i^.-i^ + l....,^-l, (B4) 
where the commutation relations: 

[Qr,Pr'] = ih6rr', [^r, = [Pr,P*r'] = 0, (B5) 

hold with q* = q^r and p* = 

By the canonical transformations given by Eqs. (|T8l) . f lTOj) . fIBip and flB2p . if/ in Eq. (jl]) 
becomes 

= V 5^ 0:4 + ^ E ^t-cV'iv^Qogo. (B6) 

s=-Afs/2 r=-7Vs/2 
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Equations of motion for Qk and are expressed by 

MQ, = -Mnl Qs for s ^ 0, (B7) 



with 



MQo = -Mnl Qo + c^NsNb % for s = 0, (B8) 
mqr = —mCjl qr for r 7^ 0, (B9) 

mgo = -mul % + c\/ NsNb Qo for r = 0, (BIO) 



MQ'i = D + 4K sin' {^] +Nbc, (Bll) 



4k sin^ [ — ]+Nsc. (B12) 

Nb, 



Equations (IB7|) - flB10P show that Qo and go are coupled although they are decoupled from 



the rest of variables. From Eqs. (1B7I) - (1B10I) . we obtain eigenfrequencies which have been 
presented in Sec. Ill B. 
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FIG. 1: Temperature dependences of energy and specific heat of a coupled HO system and an 
uncoupled HO bath with cq = 0.0. (a) Eb/Nb (dashed curve), and Eg/Ng for various Ng- 
Ng = 4 (chain curves), 10 (dotted curves), 100 (double-chain curve) and 1000 (soHd curves), (b) 
CB/ksNs (dashed curve), and Cg/ksNg for various Ng which are same as in (a), the inset showing 
Cg/ksTNg (eo = hw^i). 

FIG. 2: Temperature dependences of Cg/ksNg of coupled HO systems subjected to uncoupled 
HO baths with Ng = 4 and Nb = 100 for various cq: cq = 0.0 (solid curve), 0.01 (dotted curve), 
0.1 (dashed curve), 0.5 (bold solid curve), 1.0 (chain curve), 2.0 (bold dashed curve) and 5.0 
(double-chain curve). 

FIG. 4: Temperature dependences of Cg/kBNg of coupled HO systems subjected to uncoupled 
HO baths with Ng = 4 for various Nb'- Nb = 4 (solid curve), 10 (dotted curve), 100 (dashed 
curve) and 1000 (chain curve) with cq = 0.1, 1.0 and 5.0. Results for Nb = 4, 10, 100 and 1000 
with Co = 0.1 are indistinguishable. 

FIG. 5: Temperature dependences of C/kB (chain curve), CB/kB (dashed curve) and Cg/kB (solid 
curve) of coupled HO systems subjected to uncoupled HO baths for cq = 5.0, Ng = 4 and Nb = 10: 
C and Cb denote specific heats of system-plus-bath and bath, respectively, and the system specific 
heat is given by C5 = C — Cb- The inset shows eigenfrequencies ui for cq = 0.0 (open circles) and 
Co = 5.0 (filled circles), dashed lines being plotted for a guide of the eye (see text). 

FIG. 3: Temperature dependences of Cg/ksNg of coupled HO systems subjected to uncoupled 
HO baths for (a) co = 0.1, (b) co = 1.0 and (c) co = 5.0 with Nb = 100 for various A^^: Ng = A 
(solid curve), 10 (dashed curve) and 100 (chain curve). 
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FIG. 6: Temperature dependences of specific heat of uncoupled HO systems subjected to uncoupled 
HO baths with ujq = 1.0, D = 1.0 and K = 0.0. (a) Cs/ksNs with cq = 5.0 and Nb = 100 for 
various Ns'. Ns = 1 (solid curve), 2 (dashed curve), 4 (dotted curve) and 10 (chain curve), (b) 
C/ks (chain curve), Cs/kB (dashed curve) and Cs/kB (solid curve) for cq = 5.0, Ns = 4 and 
Nb = 10: the inset shows eigenfrequencies z/j for cq = 0.0 (open circles) and cq = 5.0 (filled circles), 
dashed lines being plotted for a guide of the eye (eo = Hujo). 



FIG. 7: Temperature dependences of specific heats of coupled HO systems subjected to coupled 
HO baths expressed by Eq. ([60|) with k = 1.0, K = 1.0 and D = 0.0. (a) Cs/ksNs with cq = 5.0 
and Nb = 100 for various Ns'- Ns = 4 (solid curve), 10 (dashed curve) and 100 (chain curve), 
(b) C/kB (chain curve), CB/kB (dashed curve) and Cs/kB (solid curve) for cq = 5.0, Ns = 4 and 
Nb = 20: the inset shows eigenfrequencies Ui for cq = 0.0 (open circles) and co = 5.0 (filled circles), 
dashed lines being plotted for a guide of the eye (eo = hyk/m). 
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